NPS  ARCHIVE 
1968 
GALLO,  S. 


COMPARISON    OF   MEASURED   AND   CALCULATED 
FLOW   PROPERTIES   AFTER  THE    ROTATING   CASCADES 
A  THREE-STAGE  AXIAL-FLOW  COMPRESSOR 


Salvatore   Frank  Gallo 


WAVAL  P0S1 

GALXFo  93 


UNITED  STATES 
NAVAL  POSTGRADUATE  SCHOOL 


THESIS 


COMPARISON  OF  MEASURED  AND  CALCULATED 

FLOW  PROPERTIES  AFTER  THE  ROTATING  CASCADES  OF  A 

THREE -STAGE  AXIAL-FLOW  COMPRESSOR 

by 

Salvatore  Frank  Gallo 


June  1968 


This  document  is  subject  to  special  export  co 
trols  and  each  transmittal  to  foreign  government 
or  foreign  nationals  may  be  made  only  with  prior, 
approval  of  the  U.    S.    Naval  Postgraduate  School. 


DIXJICYKNOXLibKARY 
NAVAL  POSTGRADUATE  SCHOOL 
MONT0t£Y  CA   93943-5101 


COMPARISON  OF  MEASURED  AND  CALCULATED 

FLOW  PROPERTIES  AFTER  THE  ROTATING  CASCADES  OF  A 

THREE-STAGE  AXIAL-FLOW  COMPRESSOR 

by 

Salvatore  Frank  Gallo 

Lieutenant  Commander,  United  States  Navy 

B.  S.,  State  University  Maritime  College,  New  York,  1957 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

AERONAUTICAL  ENGINEER 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1968 


ABSTRACT 

This  work  was  undertaken  to  check  the  feasibility  of  obtaining  the 
velocity  profiles  downstream  of  the  rotors  of  turbomachines  by  measuring 
only  the  radial  distributions  of  total  temperature,  total  pressure,  and 
absolute  flow  angle  and  using  these  data  to  calculate  the  flow  velocities 
by  the  non-isentropic  radial  equilibrium  equation  such  that  the  measured 
flow  rate  is  matched.   Curvature  effects  are  ignored  in  this  analysis. 

The  velocities  calculated  by  this  method  are  compared  with  those 
obtained  by  measurements  with  flow  probes  and  hot-wire  anemometers. 

The  tests  were  conducted  on  a  3 -stage  axial  flow  compressor  of  the 
Turbo-Propulsion  Laboratory  at  the  Naval  Postgraduate  School.  Measure- 
ments were  made  over  the  entire  operating  range  of  the  compressor  at  a 
fixed  rotor  speed  of  about  1200  rpm. 
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SECTION  1 
INTRODUCTION 

The  experimental  analysis  of  performance  of  turbomachinery  requires 
detailed  and  accurate  measurement  of  the  flow  properties  within  the  pas- 
sages of  the  machine  under  investigation. 

Even  at  the  present  state  of  the  art  it  is  necessary  to  make  a  num- 
ber of  simplifying  assumptions  for  the  design  of  turbomachines.   A  limited 
amount  of  basic  design  data  is  available  which  may  not  apply  directly  to 
a  particular  type  of  machine.   It  is  necessary,  therefore,  to  check  the 
accuracy  of  the  design  method  by  testing  actual  machines.   For  axial 
type  turbomachines  it  is  of  great  interest  to  know  the  flow  profiles 
after  the  different  stages  since  this  information  indicates  whether  the 
flow  deflections,  losses,  and  velocity  distributions  on  which  the  design 
was  based  have  been  realized.   Comparison  of  assumed  and  actual  flow, 
properties  will, make  it  possible  to  correct  or  improve  the  design  pro- 
cedure with  the  result  that  the  performance  of  later  machines  can  be 
predicted  with  greater  accuracy. 

The  measuring  of  the  flow  properties  between  the  rows  of  blades  of 
axial  turbomachines  is  connected  with  considerable  difficulties.   It  is 
necessary  to  use  small  probes  to  avoid  disturbing  the  flow;  and  some  of 
the  more  accurate  probes,  such  as  Pitot-static  probes,  cannot  be  arranged 
because  of  small  axial  blade  clearances.   Flow  angles,  total  pressures, 
and  total  temperatures  can  be  measured  with  relative  ease  and  sufficiently 
small  probes.   To  establish  the  velocity  profiles,  however,  it  is  neces- 
sary to  know  the  static  pressure  and  the  static  temperature.   The  meas- 
uring of  these  quantities  between  the  rows  of  turbomachines  can  usually 
be  carried  out  only  with  probes  that  have  been  calibrated  for  different 
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Mach  number,  Reynolds  number,  proximity  of  walls,  and  flow  angles. 
For  these  reasons  the  accuracy  of  such  measurements  is  low,  and  the 
experimental  procedure  is  quite  difficult. 

The  purpose  of  this  study  is  to  determine  whether  it  is  possible 
to  use  the  first-mentioned  measurements  together  with  certain  theoreti- 
cal flow  relations  to  obtain  the  flow  profiles  with  a  sufficient  degree 
of  accuracy.   These  conditions  are  checked  on  a  three -stage  axial  com- 
pressor of  large  size  where  the  axial  clearances  between  the  rows  of 
blades  make  it  possible  to  measure  the  velocities  directly  with  cali- 
brated flow  probes  and  hot-wire  instruments.   This  compressor  will 
henceforth  be  called  the  O.N.R.  compressor  since  the  Naval  Postgraduate 
School  received  this  machine  from  the  Office  of  Naval  Research,  Depart- 
ment of  the  Navy. 

The  author  greatfully  acknowledges  the  guidance  and  counsel  of 
Professor  M.  H.  Vavra  of  the  Department  of  Aeronautics  of  the  Naval 
Postgraduate  School.  Professor  Vavra  not  only  guided  this  research 
project  but  also  gave  the  author  an  appreciation  of  what  it  means  to  be 
an  engineer. 

The  author  also  wishes  to  express  his  appreciation  to  the  Office  of 
Naval  Research  which  made  the  compressor  used  in  this  study  available 
to  the  Naval  Postgraduate  School. 

Thanks  are  due  also  to  Mr.  J.  E.  Hammer  and  Mr.  R.  W.  Savage  of 
the  technical  staff  of  the  Turbo-Propulsion  Laboratory  as  well  as  to 
Mr.  T.  B.  Dunton  and  the  staff  of  the  Instrument  Shop  of  the  Department 
of  Aeronautics. 
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SECTION  2 
BASIS  FOR  THE  ANALYSIS 

For  any  flow  condition  there  are  two  principal  equations  of  conserva- 
tion which  must  be  satisfied,  namely,  the  equation  of  motion  and  the  equa- 
tion of  continuity. 

The  following  simplifying  assumptions  are  made  for  the  present  treat- 
ment: 

1.  The  fluid  is  assumed  to  be  a  perfect  gas. 

2.  The  general  flow  equations  are  applied  to  compute  fluid  proper- 
ties between  blade  rows  where  blade  forces  are  nonexistent. 

3.  The  flow  is  assumed  to  be  steady. 

4.  The  flow  is  assumed  to  be  axisymmetric.   By  this  assumption, 
all  partial  derivatives  with  respect  to  the  peripheral  coordinate  are 
taken  to  be  equal  to  zero.   Basically,  this  is  the  equivalent  of  saying 
that  the  three  dimensional  flow  through  a  blade  row  can  be  represented  by 
the  circumferential  average  of  the  flow  conditions  at  each  radial  posi- 
tion. 

5.  The  flow  is  assumed  to  be  adiabatic. 

6.  The  shearing  stresses  at  the  annulus  boundaries  are  ignored  in 
the  equation  of  motion,  but  their  effects  are  reflected  by  the  use  of 
measured  data. 

7.  The  specific  heat  of  the  gas  is  assumed  to  be  constant. 

With  these  simplifying  assumptions,  the  equation  of  motion  as  defined 
by  Vavra  (1)  and  the  equation  of  continuity,  respectively,  are 


VH  =  TVS   +  V  x  (V*V  ) 


(i) 
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v-ce  v)  =  o  (2) 

Since  this  study  deals  with  the  radial  variations  of  fluid  properties, 

it  suffices  to  use  the  radial  components  of  Eq.  1.   In  consistent  engi- 

/ 
nee ring  units, 

J%LPZk        Wdx      V*Tx~     *-7k     V«37"   (3> 

A  further  simplifying  assumption  is  made  inasmuch  as  the  radial 
accelerations  due  to  streamline  curvature  are  ignored,  thus  eliminating 
the  last. term  in  Eq.  3.    This  assumption  is  valid  for  lightly  loaded 
machines  where  the  curvatures  of  the  meridional  flow  path  are  small  as 
described  by  Harrison  (2)  and  Giamati  and  Finger  (3).   Since  both  condi- 
tions are  met  by  the  O.N.R.  compressor,  the  error  introduced  by  this 
assumption  will  be  small. 

In  consistent  engineering  units  the  equation  of  continuity  can  be 
expressed  by 

W    =  ZT%  J       PVa;td/T  Clbm/sec)  (4) 

AhUb 
The  integration  of  Eq.  4  within  the  boundary  layer  of  the  annulus 

walls  requires  extensive  data  in  the  boundary  layer.   These  data  were 

not  available  in  this  study  due  to  the  physical  dimensions  of  the  probes 

which  were  used  to  obtain  the  surveys.   The  introduction  of  a  gross 

weight  flow  blockage  factor  in  the  equation  simulates  the  effect  of  the 

boundary  layer  by  increasing  slightly  the  axial  velocity  components  to 

meet  the  required  continuity  condition.   With  the  weight  flow  blockage 


factor  B  ,  Eq.  4  becomes 


(5) 


/f  K11K 


hub 
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Equations  3  and  5  require  iterative  solutions  to  satisfy  both  equa- 
tions.  In  addition,  the  non-linearities  in  Eq.  3  require  numerical 
solutions;  hence,  the  use  of  a  high  speed  computer  is  almost  a  necessity 
to  carry  out  the  procedure. 

Equation  3  was  multiplied  by  dr,  and  integration  yields 


jjjCP§d.=J%T^-d^/v^ 


h.         2 
Vu 

n 


A. 


^*"*  J  V^TaA* 


(6) 


Stations  r  and.r.  are  two  closely  adjoining  radial  stations. 

Some  of  the  terms  in  Eq„  6  may  be  integrated  directly,  but  the 

2 

integrals  containing  entropy  changes  and  the  one  involving  V   /r  must 

be  determined  numerically. 

Across  any  blade  row  the  difference  in  entropy  at  a  given  stream- 


line may  be  expressed  by  ^ 


WrJflD.'"-  ju 


CWTViI 


(P^/PO 


(7) 


The  symbol  S  represents  the  stagnation  entropy.   Subscript  2  refers 
to  the  downstream  position  and  subscript  1  to  the  upstream  position. 
The  entropy  gradient  across  the  blade  row  on  two  closely  adjoining  stream- 
lines may  now  be  calculated.   For  descriptive  purposes  two  such  stream- 
lines will  be  subscripted  by  i  and  j. 

Subtracting  the  axial  entropy  change  across  the  blade  row  on  stream- 
line i  from  that  on  j  yields 


Ht-KiVilfca- 


(8) 


JL 


(Ttt,/T*,,) 


4- 


(■*./*.») 


j.r(WT.,,,) 


A, 
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The  numerical  value  of  the  entropy  gradient  in  the  radial  direc- 
tion may  now  be  simplified  by  assuming  that  position  1  is  far  upstream 
where  no  radial  entropy  gradient  exists   and  where  radial  gradients  of 

total  temperature  and  total  pressure  are  zero.   Therefore,  T    equals 

fcli 
T     and  P    equals  P    .   Then  Eo.  8  simplifies  to 
'lj'      Cli         Clj 


5*  ~  $ii  ~     J  ^ 


(\yjt,j 


(RV^l 


(9) 


The  integration  of  the  term  in  Eq.  6  that  contains  the  entropy 
gradient  ^S/c)r   is  then 
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The  static  temperature  appearing  in  Eq.  10  may  be  expressed  by 


the  stagnation  temperature  T  and  the  absolute  velocity  V,  or  by 


T  =    T*  - 

With  Eq.  11,  Eq.  10  becomes 


T-l     V 


(ii) 


hTirJ«=  2  4§^][*R(Tt+T'^(v!+v;)]  (12) 


The  non-linear  term  in  Eq.  6  containing  V  /r  can  be  evaluated 


numerically  by 
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For  simplification,  let 


F   =  U 


CR/HJ 


(14) 


Since  V  at  an  arbitrary  radius  is  to  be  determined,  V  and  V  in 
a  u 


Eq.  6  are  expressed  by  V  and  o(    ,  or 

EL 


(15) 


Equation  6  may  now  be  solved  for  V  at  an  arbitrary  radius  r, 

3 


\C  -  frrrs 


*  t^le<C2-^-F  ^-    J 


[2^R#r(T«-i;i)-F?R(T+T«/:)  + 

V^fF^  +|)^V4;^'^(F^*2-^)]       (16) 

With  measured  values  of  total  temperature,  total  pressure  and  flow 

angle  as  a  function  of  radius,  Eq„  16  could  be  solved  at  all  radial 

positions  if  the  reference  axial  velocity  V   were  known. 

a. 

i 

In  Eq«  5  the  density  is  a  function  of  velocity.   For  a  perfect  gas 


%r 


e. 


RT 


(17) 


From  isentropic  relationshi] 
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Combining  eqs.  17  and  18  with  the  static  temperature  of  Eq.  11 
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This  quantity  introduced  into  Eq.  5  establishes  the  flow  rate  for 

the  assumed  velocity  V   .   If  it  does  not  agree  with  the  measured  flow 

a. 

1 

rate,  a  new  reference  axial  velocity  V   is  assumed  to  obtain  an  improved 

a. 

1 

solution  of  Eq.  16;  and  the  process  is  repeated  until  agreement  between 
the  measured  and  calculated  mass  flow  rates  is  reached. 

Since  Eq.  16  contains  differences  instead  of  differentials,  it  is 
evident  that  small  increments  in  radius  are  necessary  to  insure  accuracy 
in  the  final  results. 

The  described  method  of  velocity  determination  is  similar  to  that 
presented  in  Ref.  3  for  the  design  and  off -design  analysis  of  axial  flow 
compressors. 
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SECTION  3 
COMPUTER  TECHNIQUE  FOR  SOLUTION 

A  computer  program  for  the  calculation  of  the  radial  flow  properties 
downstream  of  a  rotor  row  was  prepared  for  the  IBM-360-67  computer  in 
Fortran  IV  language. 

The  program  is  specifically  set  up  for  the  0,N.R.  three-stage  axial 
compressor  of  the  Naval  Postgraduate  School  described  in  Section  4,  but 
slight  modifications  make  it  applicable  to  radial  surveys  downstream  of 
arbitrary  axial  rotor  rows  of  either  turbines  or  compressors  for  the 
assumptions  stated  in  Section  2. 

3.1     Non-Isentropic  Radial  Equilibrium  Equation 

The  computer  program  has  been  set  up  for  data  from  a  maximum  of 
fifty  radial  stations,  not  necessarily  extending  to  the  hub  and  to  the 
tip  radii,  although  accuracy  would  be  reduced  if  data  were  not  measured 
close  to  the  hub  and  the  tip  of  the  blading.   Values  of  total  pressure, 
total  temperature,  and  absolute  flow  angle  are  introduced  as  a  function 
of  radius.   Other  input  values  are  measured  flow  rate,  initial  reference 
value  of  axial  velocity,  assumed  blockage  factor,  desired  accuracy  in  flow 
rate  matching,  and  rotor  speed. 

The  introduced  data  which  are  measured  at  discrete  radii  are  first 
expressed  with  fourth  order  polynomials  as  a  function  of  radius.   The 
polynomial  fit  for  the  data  was  obtained  by  the  method  of  least  squares. 
Due  to  the  smooth  change  of  the  data,  the  polynomial  fit  was  extremely 
accurate  throughout  the  radial  span  for  all  the  parameters.  An  example 
of  the  accuracy  of  the  curve  fit  is  given  in  Table  B-2  where  the  symbol 
X  signifies  radius;  F2,  the  input  data  for  different  radii;  Y,  the  values 
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of  the  parameter  obtained  from  the  polynomial;  and  DELY,  the  difference 
between  the  measured  data  and  those  from  the  polynomial.   The  subprogram 
used  for  the  polynomial  curve  fitting  is  available  in  the  Naval  Post- 
graduate School  subroutine  library  and  is  presented  in  Appendix  B. 

A  maximum  of  one  hundred  points  may  be  specified  at  which  the 
equilibrium  equation  is  to  be  solved.   For  brevity,  the  non-isentropic 
radial  equilibrium  equation  shall  be  referred  to  as  the  NISRE  equation 
in  the  remainder  of  this  presentation.   The  radial  positions  of  the 
points  are  obtained  by  dividing  the  radius  from  hub  to  tip  into  equal 
increments.   Values  of  total  pressure,  total  temperature  and  absolute 
angle  of  flow  are  then  obtained  from  the  polynomials  at  these  specified 
radial  positions.   The  measured  mass  flow  rate  is  calculated  from  meas- 
ured conditions  in  the  inlet  to  the  compressor. 

The  subprogram  dealing  with  the  NISRE  equation  then  solves  for  the 
axial  velocity  and  density  profile  from  blade  hub  to  blade  tip.   The 
profiles  thus  obtained  are  used  in  the  continuity  equation  which  is 
integrated  by  the  trapezoidal  method.   The  initial  axial  velocity  at 
the  reference  position  is  then  changed  arbitrarily,  and  a  second  solution 
is  generated.   The  third  and  subsequent  values  of  the  reference  axial 
velocity  are  calculated  with  the  Newton-Raphson  iteration  technique. 
This  process  continues  until  the  desired  accuracy  in  mass  flow  rate  is 
obtained. 

With  the  axial  velocity  profile  thus  determined,  all  remaining  fluid 
properties  may  be  calculated.   These  properties  are  printed  out  in  dimen- 
sional and  non-dimensional  form.   Examples  of  the  print  out  of  the  solu- 
tions are  presented  in  Tables  B-3  and  B-4. 
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3.2  Determination  of  Velocity  Profiles  with  Flow  Probes 

The  measured  values  of  total  pressure,  static  pressure,  absolute 
flow  angle,  total  temperature  and  rotor  speed  are  used  to  obtain  the 
axial  velocity  profile.   The  fluid  properties  at  different  radial 
stations  are  calculated  and  printed  out  in  a  format  similar  to  that  of 
the  NISRE  equation  (See  Tables  B-5  and  B-6.). 

3.3  Determination  of  Velocity  Profiles  with  the  Hot-Wire  Anemometer 
The  hot-wire  anemometer  survey  was  made  at  the  same  positions  as 

the  flow  probe  and  total  temperature  surveys.   To  place  the  hot-wire  as 
nearly  perpendicular  as  possible  to  the  flow,  the  absolute  flow  angle 
at  each  radial  position  was  obtained  from  a  five-hole  probe.   The  hot- 
wire output  was  set  to  a  reference  value  at  the  mean  square  radius  of 
the  annulus,  and  all  values  were  measured  with  respect  to  this  setting. 
The  velocity  determined  by  the  flow  probe  at  the  mean  square  radius  was 
used  as  the  reference  velocity  for  the  hot-wire  anemometer.   The  root 
mean  square  radius  was  used  as  reference  radius  since  the  wall  inter- 
ference effects  on  the  flow  probe  are  small  at  this  location. 

Tables  B-7  and  B-8  give  sample  outputs  of  the  flow  properties  at 
different  radii  calculated  from  the  measured  data.   It  should  be  noted 
that  the  static  pressures  and  static  temperatures  are  the  values  obtained 
from  the  analysis.   The  measured  static  pressures  are  not  used  in  these 
calculations. 
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SECTION  4 
O.N.R.  THREE -STAGE  AXIAL-FLOW  COMPRESSOR 

The  compressor  installation  is  shown  in  Fig.  1.   A  detailed  descrip- 
tion of  the  test  installation  is  given  by  Bowen  £t  al_.  (4).  Marshall 
(5)  describes  the  test  set  up  at  the  Naval  Postgraduate  School. 

The  compressor  is  suitable  for  a  large  variety  of  research  experi- 
ments.  Each  one  of  the  rotor  and  stator  blade  rows  may  be  removed. 
Although  stagger  angles  are  adjustable ,  they  were  not  varied  during  the 
present  experiments.   All  three  rotor  blade  cascades  were  set  at  a 
stagger  angle  of  43.8  degrees,  with  all  stator  blade  stagger  angles  at 
39.8  degrees. 

The  casing  of  the  compressor  was  designed  to  give  maximum  accessi- 
bility for  measuring  instruments.   Six  rectangular  instrument  ports  are 
located  in  the  upper  half  of  the  casing.   These  ports  accommodate  a 
special  instrument  carriage  permitting  detailed  surveys  in  radial  and 
peripheral  directions  (Fig.  2).   Three  of  the  instrument  ports  are 
suitable  for  stator  outlet  surveys  over  a  peripheral  arc  of  15  degrees, 
and  the  other  three  are  suitable  for  rotor  outlet   surveys  over  the 
same  arc.   Only  rotor  outlet  surveys  were  made  during  this  study. 

The  casing  has  numerous  radial  survey  ports.   Several  of  these 
were  used  downstream  of  each  rotor. 

The  compressor  is  electrically  driven  and  operates  at  two  fixed 
speeds  of  either  about  900  rpm  or  1200  rpm. 

The  inlet  duct  (Fig.  3)  consists  of  a  screen,  an  entrance  bellmouth 
and  a  length  of  straight  pipe  with  a  honeycomb  flow  straightener  ahead 
of  the  inlet  guide  vanes. 
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The  compressor  discharge  (upper  foreground  of  Fig.  1)  is  fitted 
with  an  electrically  operated  throttle  valve  whose  positions  are 
indicated  by  a  counter. 
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SECTION  5 


INSTRUMENTATION 


The  compressor  inlet  duct  allows  vertical  and  horizontal  surveys 
to  be  made  across  the  inlet  as  well  as  at  60  degrees  from  the  vertical 
direction  by  unbolting  and  rotating  the  inlet  duct.  A  3/32  inch  diam- 
eter Prandtl-type  Pitot-static  probe  was  mounted  for  vertical  survey 
(Fig.  4)  with  readout  on  a  Meriam  Instrument  Company  Model-750  Micro- 
manometer  (Fig.  5).  A  1/4  inch  diameter  hot-wire  anemometer  probe  was 
mounted  for  horizontal  survey.   The  hot-wire  probe  was  constructed  at 
the  Naval  Postgraduate  School  and  provides  for  quick  changing  of  tips 
in  the  event  of  wire  breakage  (Fig.  6).  A  supply  of  tips  was  kept 
readily  available  during  the  inlet  surveys,  and  tip  replacement  was 
effected  on  several  occasions  while  the  compressor  was  in  operation. 
Hot-wire  data  were  read  on  a  Security  Associates  Applied  Sciences  Model 
200  Constant  Temperature  Anemometer.   Details  of  the  flow  rate  measure- 
ment are  given  in  Section  6. 

Wet  and  dry  bulb  Mercury  glass  thermometers  were  mounted  on  the 
inlet  screen. 

A  1/4  inch  diameter  Prandtl-type  Pitot-static  probe  was  installed 
between  two  of  the  inlet  supports  ahead  of  the  guide  vanes  (foreground 
of  Fig.  7)  with  readout  on  the  same  micromanometer  as  the  inlet  duct 
Pitot-static  probe.   This  was  made  possible  by  the  use  of  two  scanner 
valves  (Fig.  8).   In  addition,  a  1/8  inch  diameter  Kiel  probe  with  a 
total  type  thermocouple  was  mounted  in  the  same  meridional  plane  through 
one  of  the  radial  access  ports.   This  total  temperature  probe  is  shown 
in  the  foreground  of  Fig.  7.   The  reference  junction  of  the  thermocouple 
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was  placed  in  an  ice  bath.   The  output  of  the  thermocouple  was  read  on 
a  Leeds  and  Northrup  Potentiometer  (Fig.  9). 

Radial  surveys  were  conducted  downstream  of  each  of  the  three  rotor 
blade  rows.  Total  and  static  pressure  measurements  were  obtained  with 
a  United  Sensor  Corporation  DA120  five  hole  probe  (serial  number  535) 
although  the  pitch  angle  was  assumed  to  be  zero.   The  five  hole  probe 
was  mounted  in  the  special  carriage  mentioned  in  Section  4  with  readout 
on  a  water  manometer  board  (Fig.  10) .   Flow  direction  was  obtained  by 
balancing  pressure  on  the  static  ports.   Total  pressure  and  the  differ- 
ence between  total  pressure  and  static,  pressure  were  measured  with  the 
same  manometer  board.  Total  temperature  measurements  at  different  radii 
were  obtained  by  use  of  a  thermocouple  mounted  in  a  1/4  inch  diameter 
Kiel  probe  which  was  placed  in  one  of  the  radial  access  holes.   The 
reference  junction  for  this  thermocouple  was  also  placed  in  an  ice  bath, 
and  the  aforementioned  potentiometer  was  used  for  the  readout  by  the  use 
of  the  double-pole  double-throw  knife  switch  shown  in  the  background  of 
Fig.  9.   Hot-wire  anemometer  data  at  different  radii  were  measured  after 
each  of  the  three  rotor  blade  rows  by  means  of  a  1/4  inch  diameter  hot- 
wire probe.  Figure  11  shows  the  installation  of  the  total  temperature 
probe  and  the  hot-wire  anemometer  probe  used  for  the  radial  surveys.   The 
hot-wire  anemometer  readout  instrument  is  shown  in  Fig.  9. 

Compressor  speed  was  measured  with  a  flux  cutter  mounted  on  the 
compressor  drive  shaft  whose  impulses  were  read  on  a  Monsanto  Model  100A 
Counter- timer. 
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SECTION  6 
FLOW  RATE  CALIBRATION 

The  purpose  of  the  flow  calibrations  was  to  obtain  a  factor  which, 
when  multiplied  by  the  velocity  ahead  of  the  inlet  guide  vanes,  gives 
the  volume  flow  rate  through  the  compressor, 

Marshall  (5)  found  a  flow  rate  calibration  factor  for  the  O.N.R. 
compressor  by  means  of  horizontal  and  vertical  surveys  of  the  inlet  with 
the  existing  Pitot-static  probe.  This  calibration  factor  was  to  be 
checked  by  a  more  accurate  survey  with  a  hot-wire  anemometer. 

Since  Marshall  showed  that  the  velocity  profiles  in  the  vertical 
and  horizontal  planes  were  similar,  it  was  decided  to  conduct  a  survey 
in  the  horizontal  plane  only. 

Since  hot-wire  instruments  give  only  the  ratio  of  the  actual  velocity 
to  some  reference  velocity,  initial  hot-wire  instrument  setting  was  made 
at  the  inlet  duct  centerline  where  the  velocity  obtained  from  readings  of 
the  Prandtl  type  Pitot-static  probe  was  used  as  the  reference  velocity. 
At  the  inlet  duct  centerline  the  Pitot-static  probe  is  18  inches  from  the 
duct  wall;  hence,  wall  effects  are  minimized.   The  hot-wire  instrument 
was  set  at  90  per  cent  of  full  readout  range  at  the  centerline,  and 
surveys  were  made  at  radial  increments  of  3  inches  to  a  radius  of  15 
inches  and  of  0.2  inches  to  a  point  0,2  inches  from  the  duct  wall.   The 
design  of  the  hot-wire  probe  prevented  measurements  closer  to  the  wall. 
Frequent  checks  were  made  at  the  duct  centerline  to  ascertain  that  the 
centerline  setting  of  90  per  cent  had  not  drifted.   Such  instrument 
drifts  were  noticed  only  on  rare  occasions .   If  drift  occurred,  it  was 
always  less  than  2  per  cent  of  full  scale,  and  the  hot-wire  instrument 


30 


was  then  i-eset  to  90  per  cent.   Surveys  of  this  type  were  conducted  on 

both  sides  of  the  centerline.   Total  pressure  and  velocity  head  readings 

were  also  taken  at  the  duct  centerline  and  ahead  of  the  inlet  guide 

vanes  with  the  fixed  Pitot-static  pressure  probe  connected  to  a  water 

micromanometer.  Wet  and  dry  bulb  temperatures  were  recorded  from  the 

Mercury  thermometers  on  the  inlet  screen.   The  temperature  in  the  test 

cell  was  recorded  also  for  the  purpose  of  correcting  water  column 

heights . 

The  following  symbols  are  used  for  the  recorded  data: 

qlv       -velocity  head  at  the  centerline  of  the  inlet  duct  (in. HO) 

P_       -total  pressure  at  the  centerline  of  the  inlet  duct  (in.H  0) 
t,  2 

lr 

q~       -velocity  head  ahead  of  the  inlet  guide  vanes  (in.  HO) 

P        -total  pressure  ahead  of  the  inlet  guide  vanes  (in.  HO) 

d2r  l 

TOR      -temperature  at  the  inlet  screen  (°F) 

BAR       -barometric  pressure  (in.Hg.) 

TBAR      -temperature  of  the  Mercury  column  (°F) 

TW       -temperature  of  the  water  column  (*F) 

VRR      -hot-wire  anemometer  readings,  right  of  centerline  looking 

downstream  (dimenslonless) 
VRL      -hot-wire  anemometer  reading,  left  of  centerline  looking 

downstream  (dimenslonless) 
The  specific  gravity  of  the  Mercury  in  the  barometer  was  corrected 
for  temperature  by 

GHGRM  =  l3.63305"-0.00l3J.303O3TBAR    (20) 
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The  specific  gravity  of  the  water  column  was  corrected  for  tempera- 
ture by 

GWR [1=039837633+1. 0605" r5£TW*IO~4-!.S"9 31 86TW**  J0~6     <21> 

These  specific  gravity  relations  were  obtained  from  tabulated  data 
of  Ref.  6. 

With  these  corrections,  pressures  were  converted  to  absolute  values, 
or 

PAKIB=  BAR  C0.4S9I5"85")  GHGRM  j^-     (  Psf<0  (22) 

where  the  constant  0, 4891585  is  the  conversion  factor  from  inches  of 

Hg.  to  pounds  per  square  inch  for  a  specific  gravity  of  Mercury  of  13.54. 

Then 

pbx  =  PtJ  G-WRM)  ^l13Z    +  PAMB         (psf<^)  (23) 

The  constant  62.42732  is  the  conversion  factor  from  feet  of  water 
to  pounds  per  square  foot  for  a  specific  gravity  of  water  of  1.0. 


The  static  pressure  P  was  then  calculated  by 

P.  =  p*-#,  (psfa.) 


(24) 


(25) 


Because  of  the  low  velocities  of  about  70  feet  per  second  in  the 
inlet  duct,  and  ahead  of  the  compressor  inlet  guide  vanes,  compressi- 
bility effects  were  ignored,  and  the  temperature  TOR  at  the  inlet  screen 
was  used  as  the  static  temperature „ 

T.  =  TOR  +  453.7      CR)  «» 


32 


Considering  the  air  entering  the  compressor  as  a  perfect  gas  having 
a  gas  constant  R  of  53.35  (ft-lb)/(lbm-°R) ,  the  mass  density  is 

For  each  run  the  density  was  taken  as  the  average  obtained  from  all 
readings. 

The  velocity  at  the  inlet  duct  centerline  was  then  calculated  by 


V,  -   /-^  (ft/j.«)  (28) 

This  velocity  was  used  as  the  reference  velocity  for  the  hot-wire 

anemometer  setting.   Then,  the  average  velocity  V  at  an  arbitrary  radius 
is  taken  as 


v  =  (VRL    VRR)   SL     (ft/sec) 


(29) 


The  constant  0.9  must  be  introduced  since  the  hot-wire  anemometer 
was  set  at  90  per  cent  of  full  scale  readout  of  the  instrument  at  the 
centerline  in  case  velocities  higher  than  those  at  the  centerline  would 
occur. 

The  volume  flow  rate  was  calculated  by 


V0LFL=  3.14-IC  V,  (-1)  +  3.I4U  Z.W,  (n^/t^)  - 

i 


The  subscript  1  in  the  summation  refers  to  the  inlet  duct  center- 
line.   The  subscript  n  pertains  to  the  reading  0.2  inches  from  the  duct 
wall,  and  the  subscript  n+1  is  for  the  conditions  at  the  duct  wall. 
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The  velocity  V  ahead  of  the  inlet  guide  vanes  was  obtained  in  a 
similar  manner,  using  the  temperature  measured  at  the  inlet  screen  as 
the  static  temperature  since  the  difference  between  total  and  static 
temperature  is  less  than  Q„4°F  at  the  highest  velocity  of  70  feet  per 
second. 

The  calibration  constant  for  each  run  was  obtained  from 

CCC  =  -VQWL F L  (ftVsec   per  ft/sec)  OD 

3 

This  constant  was  found  to  be  4 „ 4050  ft  /sec  per  ft/sec  velocity 

obtained  ahead  of  the  inlet  guide  vanes.   The  maximum  relative  devia- 
tion of  CCC  for  all  surveys  was  0,3  per  cent. 

This  data  reduction  was  carried  out  on  the  IBM-360-67  computer.   The 
program  is  explained  in  Appendix  A  which  also  gives  a  sample  of  the 
computer  output. 
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SECTION  I 


RADIAL  SURVEYS 


Radial  surveys  were  conducted  downstream  of  each  row  of  rotor  blades 
of  the  O.N.R.  compressor.   Radial  distributions  of  total  pressure,  veloc- 
ity head,  total  temperature,  absolute  flow  angle,  and  hot-wire  velocity 
ratio  were  obtained  at  various  flow  rates. 

Flow  rate  was  varied  by  means  of  the  discharge  throttle  valve.   After 
setting  the  throttle  valve  to  a  new  position,  data  were  taken  only  after 
waiting  for  about  20  minutes  to  stabili?;e  the  temperature  in  the  test 
cell.   The  test  cell  has  a  front  door  through  which  air  is  normally  admit- 
ted to  the  compressor  from  the  atmosphere.   The  air  is  discharged  back 
into  the  atmosphere  through  a  door  located  opposite  to  the  front  door  and 
aft  of  the  compressor.  However,  experience  showed  that  gusts  in  the 
atmosphere  produced  erratic  measurements  particularly  with  the  hot-wire 
instrumentation.  By  closing  the  front  door  it  was  found  that  the  flow 
oscillations  could  be  minimised,  but  the  above-mentioned  waiting  period 
was  necessary  for  the  air  to  reach  a  new  equilibrium  temperature  in  the 
te3t  cell.   Steady  temperature  conditions  were  checked  by  the  mercury 
thermometer  at  the  inlet  screen  and  the  total  temperature  Kiel  probe 
ahead  of  the  inlet  guide  vanes. 

Periodic  readings  of  the  total  pressure,  velocity  head  and  total 
temperature  ahead  of  the  inlet  guide  vanes  were  recorded  during  each  run 
for  the  determination  of  the  flow  rates.  Conditions  at  the  inlet  screen 
and  ahead  of  the  inlet  guide  vanes  were  found  to  be  very  stable  for  all 
flow  rates  irrespective  of  throttle  settings  and  atmospheric  conditions. 

The  previously  described  traversing  instrument!:  for  the  radial 
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surveys  were  mounted  downstream  of  each  one  of  the  rotors.   Initially, 
a  peripheral  survey  was  made  with  the  flow  probe  at  various  radii  to 
check  whether  the  flow  conditions  were  uniform  in  the  peripheral  direc- 
tion.  Generally,  this  was  found  to  be  true  downstream  of  each  row  of 
rotor  blades  at  flow  rates  above  full  stall. 

For  most  of  the  radial  survey  runs  a  hot-wire  instrument  setting  of 
90  per  cent  was  used  at  a  reference  radius  of  14.5  inches,  or  very  close 
to  the  root  mean  square  radius  of  the  annulus.   To  properly  position  the 
hot-wire  with  respect  to  flow  direction,  initial  angle  setting  of  the 
probe  was  made  with  the  flow  angles  obtained  from  the  flow  probe.   This 
setting  was  then  refined  by  turning  the  hot-wire  into  the  direction 
which  gave  maximum  output.  Almost  without  exception,  the  hot-wire  angle 
position  thus  found  was  identical  with  the  angle  measured  by  the  flow 
probe.   The  instrument  setting  of  90  per  cent  at  the  radius  of  14.5 
inches  was  periodically  checked  for  drift.   Drift  was  encountered  only 
occasionally  and  was  always  less  than  2  per  cent.   If  drift  was  noticed, 
the  hot-wire  instrument  was  reset  to  90  per  cent  of  full  scale  at  the 
reference  radius. 

The  total  temperature  Kiel  probe  was  set  at  the  angle  which  was 
obtained  from  the  flow  probe  at  the  particular  radius. 

For  each  run,  radial  survey  readings  of  total  temperature,  total 
pressure,  velocity  head,  absolute  flow  angle,  and  hot-wire  anemometer 
velocity  ratio  were  taken  from  the  blade  tip  radius  of  18  inches  to  a 
distance  of  0.2  inches  from  the  hub  which  has  a  radius  of  10.8  inches. 
Although  readings  closer  to  the  hub  would  have  been  desirable,  the  design 
of  the  flow  probe  prevented  readings  at  distances  less  than  0.2  inches 
from  the  hub. 
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Data  were  taken  at  radial  increments  of  0.5  inches  for  the  majority 
of  the  runs.   Initially,  the  increment  in  radius  was  chosen  as  0.2  inches. 
However,  because  of  the  smooth  variation  of  the  flow  conditions  in  the 
annulus,  this  increment  could  be  increased  to  0.5  inches  without  loss 
of  accuracy. 

Various  flow  rates  were  set  by  means  of  the  discharge  throttle  valve. 
The  throttle  valve  position  is  indicated  by  means  of  a  digital  counter. 
The  counter  is  calibrated  from  0  to  1360  (closed  to  open).   Radial  surveys 
downstream  of  each  rotor  row  were  taken  at  counter  readings  of  1000,  900, 
800,  700,  600,  500,  and  435.   Below  a  counter  reading  of  435,  low  rumbling 
noises  and  vibrations  were  observed  as  well  as  significant  read-out  fluc- 
tuations, indicating  that  major  portions  of  the  blading  were  stalled. 
Peripheral  surveys  with  the  flow  probe  at  a  counter  setting  below  435 
showed  significant  flow  variations  with  respect  to  the  peripheral  direc- 
tion.  These  conditions  existed  between  the  hub  and  approximately  mid- 
blade  height.   In  addition,  extx-eme  variation  in  the  hot-wire  anemometer 
reading  was  noticed  also  in  this  region.   In  some  cases,  the  hot-wire 
oscillations  approached  30  per  cent  of  the  instrument  read-out.   Because 
of  this  condition  no  surveys  were  made  at  throttle  settings  below  435 
counts.   At  a  throttle  setting  of  435  counts,  slight  fluctuations  in  the 
hot-wire  anemometer  readings  were  noticed  between  the  hub  radius  of  10.8 
inches  and  a  radius  of  12  inches,  whereas  the  flow  probe  readings  were 
steady  in  this  portion  of  the  annulus. 

No  flow  profiles  are  presented  for  readings  of  the  throttle  valve 
counter  above  1000,  although  four  such  surveys  were  undertaken.   They  do 
not  show  appreciable  differences  since  the  flow  rate  increases  only  by 
about  2  per  cent  for  counter  readings  between  1000  and  1360. 
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SECTION  8 
DISCUSSION  OF  RESULTS 

8.1  General 

For  the  present  study  the  O.N.R.  compressor  was  run  for  68  hours, 
of  which  28  hours  were  required  for  the  flow  rate  calibrations. 

The  tests  were  conducted  with  the  so-called  free-vortex  blading 
of  the  compressor  installed.   The  rotor  blades  of  all  three  rotor  rows 
were  set  at  the  design  stagger  angles  of  43.8  degrees  with  respect  to 
the  axial  direction.   All  stator  blades  had  stagger  angles  of  39.8 
degrees.   The  runs  were  made  at  the  high  rotor  speed  of  about  1190  rpm. 

8.2  Flow  Rate  Calibration 

The  flow  rate  calibrations  were  undertaken  to  verify  the  calibra- 
tion factor  found  by  Marshall  (5).   The  fluctuating  readings  due  to 
atmospheric  conditions  experienced  by  Marshall  were  eliminated  by  clos- 
ing the  door  in  front  of  the  compressor  inlet  duct.   It  is  interesting 
to  note  that  during  the  radial  surveys  downstream  of  the  third  rotor, 
the  door  to  the  atmosphere  was  opened.   Even  though  the  atmospheric 
conditions  were  calm,  fluctuations  of  up  to  10  per  cent  were  noticed 
in  the  hot-wire  readings.   The  fluctuations  stopped  when  the  door  was 
closed. 

Velocity  measurements  in  the  inlet  duct  were  made  from  the  center- 
line  to  a  radius  of  17.8  inches.   It  was  not  possible  to  measure  the 
velocities  closer  than  0.2  inches  from  the  wall.   Since  the  inlet  pipe 
has  a  diameter  of  36  inches,  the  hot-wire  probe  had  to  be  supported 
from  both  walls  to  eliminate  vibration.   Therefore,  the  probe  had  to 
be  attached  at  right  angles  to  the  support  structure  (Fig.  6)..   The 
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probe  thickness  of  1/4  inch  and  the  fillet  at  the  attachment  point 
prevented  measurements  closer  than  0.2  inches  from  the  wall.   Velocity 
measurements  to  the  wall  would  have  been  desirable  for  a  more  accurate 
calculation  of  volume  flow  rate.   In  the  calculation  of  volume  flow 
rate,  the  velocity  measured  at  a  radius  of  17.8  inches  was  assumed  to  be 
that  at  the  wall  also. 

8.3     Presentation  of  Data 

The  axial  velocity  profiles  obtained  from  the  NISRE  solution  and 
from  the  flow  probe  and  hot-wire  measurements  are  presented  in  Figs. 
12  through  32  in  dimensionless  form.   These  quantities  were  made  dimen- 
sionless  by  dividing  the  velocities  by  the  rotor  tip  speed.   The  parame- 
ter  CD   thus  obtained  is  convenient  for  two  reasons.   First,  assuming 

that  a  mean  axial  velocity  V  exists  in  the  flow  annulus,  the  volume 

a 

flow  rate  can  be  expressed  by 

VP=  V.,   |44      (ft-Vsec)  (32) 

and  the  rotor  tip  speed  is 

UT  -  OJ  -jf-  (ft/sec)  (33) 

where  00  is  the  angular  velocity  of  the  rotor  in  radians  per  second. 
Thus 


*    =    ffHfiQft}-  .#CW4.7« 


(34) 


Therefore,   Cp   is  a  flow  coefficient  relating  volume  flow  rate  to 
rotor  speed.   A  numerical  average  of  the  values  of   (P  obtained  from 
the  radial  distribution  of  axial  velocities  at  a  given  flow  rate  gives 
an  average  value  of  (Q      which  is  a  measure  of  the  volume  flow  rate  at 
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a  particular  throttle  setting  and  rotor  speed.   Second,  the  flow  coef- 
ficient is  a  convenient  parameter  for  comparing  the  velocities  obtained 
by  the  NISRE  solution  and  those  measured  with  the  flow  probe  and  the 
hot-wire  anemometer  because  changes  in  rotor  speed  have  no  influence  on 

The  total  pressures  measured  with  the  flow  probe  downstream  of 
each  of  the  rotors  at  various  flow  rates  are  presented  in  Figs  33  through 
35.   These  pressures  are  shown  in  dimensionless  form  by  the  parameter 

where  the  subscript  "a"  refers  to  the  atmospheric  conditions. 

The  work  input  of  a  blading  element  in  an  axial-flow  compressor 
can  be  expressed  by  the  so-called  Euler  turbine  equation.   The  rise  in 
total  temperature  from  station  (1)  ahead  of  to  station  (2)  after  the 
blading  is 

ATtr-^prUtVu.-A.Vu,)      rR)         (36) 

For  a  frictionless  process  of  an  incompressible  fluid,  the  theoretical 
rise  in  total  pressure  (  /^  P  )*  then  becomes 

(AP*)*=  pCO  OTaVu,-*,^,)     (psf)  (37) 

With  losses,  the  pressure  rise   /^  P  is 

APr=  1   (&P*f  =  ^(^Vu.-K.VJ  (p«f)  (38> 

where   ^   is  the  stage  efficiency.   Rewritten 
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The  so-called  pressure  coefficient  "Uj  is  therefore  a  function  of  the 
flow  angles  and  the  losses  only  and  is  independent  of  the  density  of 
the  flow  and  the  rotor  speed,  as  long  as  the  influences  of  Reynolds  number 
and  Mach  number  can  be  ignored.   The  aforementioned  parameter  iX1  is  suit- 
able for  comparing  the  radial  distributions  of  the  total  and  static  pres- 
sures measured  after  the  rotors  because  the  influences  of  varying  densi- 
ties and  rotor  speeds  are  eliminated. 

Since  the  present  study  was  not  concerned  with  determining  blade 
losses,  no  measurements  of  total  pressure  were  made  upstream  of  the  rotors 
while  the  radial  surveys  were  conducted  downstream  of  them.   Since  the 
surveys  after  the  rotors  were  conducted  on  different  days,  it  was  not 
possible  to  maintain  exactly  the  same  value  of  the  average  flow  coef- 
ficient  v-P  ;  hence,  the  values  of  ^IjJ   presented  for  the  total  pressure 
cannot  be  used  to  establish  the  quantity  VJ/  defined  earlier. 

The  static  pressures  calculated  by  the  NISRE  method,  those  measured 
by  the  flow  probe,  and  the  ones  calculated  by  using  the  velocity  measured 
by  ttye  hot-wire  anemometer  were  made  dimensionless  by  the  same  method. 
The  dimensionless  parameter  thus  obtained  is  convenient  for  comparing  the 
static  pressures  determined  by  these  different  methods  because  changes  of 
atmospheric  pressure,  density,  and  tip  speed  do  not  influence  the  data. 

Figures  36  through  56  show  the  measured  values  of  the  absolute  flow 
angles  0(  obtained  for  all  runs. 

8.4     Accuracy  of  the  NISRE  Method 

Figures  12  through  32  show  very  good  agreement  between  the  velocity 
profiles  calculated  by  the  NISRE  method  and  those  measured  by  the  flow 
probe  and  hot-wire  anemometer.   The  agreement  is  good  downstream  of  all 
rotors  through  the  range  of  flow  rates,  with  the  possible  exception  of 
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Fig.  24.   This  figure  shows  axial  velocity  profiles  obtained  after  the 
second  stage  rotor  at  a  flow  rate  of  22.85  lbm/sec  corresponding  to  a 
discharge  throttle  valve  counter  setting  of  500.   This  run  was  repeated 
to  check  the  results;  but  for  both  runs  they  were  almost  identical. 
The  calculated  velocities  differ  from  the  measured  ones  by  as  much  as 
7.5  per  cent  at  radius  ratios  R/Rn  of  0.67  and  0.9.   No  obvious  explana- 
tion could  be  found  for  this  deviation  from  the  otherwise  close  agree- 
ment between  the  calculated  and  measured  data. 

Figures  57  through  77  show  the  comparison  of  the  static  pressures 
obtained  from  the  three  methods.   This  is  not  an  independent  comparison 
since  these  pressures  are  dependent  on  the  velocities  also.   With  the 
flow  probe,  the  static  pressures  were  measured  directly,  but  for  the 
NISRE  calculation  and  hot-wire  measurement,  the  static  pressures  were 
calculated  by  using  Eqs.  11,  17,  and  19.   However,  good  agreement  has 
been  obtained  with  the  three  methods. 

8.5     Blockage  Factor 

The  value  of  the  gross  weight  flow  blockage  factor  to  be  used  with 
the  NISRE  method  has  a  large  influence  on  the  results.   Reference  2 
suggests  a  blockage  factor  of  0.98  for  the  first  stage  and  0.96  for  all 
subsequent  stages. 

Blockage  factors  of  0.96  for  the  second  and  third  stage  rotors  gave 
the  results  of  Figs.  19  through  32.   Blockage  factors  of  0.98  were  first 
tried  for  the  calculations  of  the  velocity  profiles  for  the  first  stage 
rotor  but  gave  velocities  well  below  those  measured  by  the  flow  probe  and 
the  hot-wire  anemometer.   The  effect  of  the  area  blockage  by  the  wall 
boundary  layers  was  then  increased  by  changing  the  blockage  factor  to  0.97 


42 


and  finally  to  0.96.   This  increase  in  area  blockage  did  not  materially 
change  the  shape  of  the  calculated  profiles  but  increased  the  magnitudes 
of  the  velocities.   Blockage  factors  of  0.96  for  the  first  rotor  gave  the 
best  correlation  between  the  axial  velocities  calculated  by  the  NISRE 
method  and  those  measured.   Figures  12  through  18  show  the  axial  velocity 
profiles  after  the  first  rotor  for  this  value.   It  is  believed  that  the 
requirement  for  a  smaller  blockage  factor  after  the  first  rotor  is  caused 
by  the  boundary  layer  build-up  along  the  walls  of  the  inlet  duct  which 
has  a  length  of  six  feet. 

8.6     Compressor  Design  Point 

As  stated  in  subsection  8.1,  the  tests  were  conducted  with  the  so- 
called  free-vortex  blading  of  the  compressor.  At  their  design  point,  such 
bladings  are  supposed  to  produce  flow  patterns  where  the  product  of  the 
radius  r  and  the  peripheral  component  V  of  the  absolute  velocity  remain 
constant  in  the  radial  direction  after  each  row  of  blades.   If  entropy 
gradients  are  ignored,  free-vortex  flows  have  constant  total  pressure  P  , 

constant  total  temperature  T  ,  and  constant  axial  velocity  components  V 

t  a 

between  the  rows  of  blades.   However,  if  entropy  gradients  occur  even 
though  total  temperatures  do  not  change  in  the  radial  direction,  the  axial 
velocity  components  of  the  absolute  velocity  are  not  constant.   Equation 
3,  without  the  term  involving  meridional  streamline  curvature,  can  be 
expressed  also  by 

J$CP  ITT  "   ¥^j   n   "37T~   +  ^  d* 
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Hence,  if  the  product  rV  were  constant  in  the  radial  direction,  the 
term  containing  V  would  be  zero.   If  the  total  temperature  were  also 
constant  in  the  radial  direction,  the  only  terms  remaining  are  those 
involving  entropy  and  axial  velocity  component.   This  indicates  that 
with  the  aforementioned  assumptions,  the  variations  of  the  axial  velocity 
components  are  a  function  of  the  entropy  gradients  and  the  static  tempera- 
tures. 

From  the  presented  data  it  can  be  noticed  that  neither  constant 

velocities  V  nor  constant  total  pressures  P  were  reached  with  the 
a  t 

installed  blading,  irrespective  of  flow  rate.   In  Ref.  3,  the  design 

stagger  angle  of  the  stator  blades  is  given  as  29.8  degrees.   Hence, 

since  a  value  of  39.8  degrees  was  used  for  the  present  tests,  it  is 

likely  that  radial  changes  in  V  and  P  were  due  to  this  condition. 

■  a      t 

All  the  axial  velocity  profiles  calculated  by  the  NISRE  method  show 
regions  of  increased  axial  velocities  at  radius  ratios  R/Rn  of  about 
0.67  and  0.90.   At  which  of  these  two  stations  the  maximum  velocity 
occurs  depends  on  the  flow  rate.   As  flow  rate  decreases,  the  region  of 
highest  velocity  moves  from  the  inner  to  the  outer  part  of  the  annulus. 
The  measured  velocity  profiles  show  the  same  trend  but  are  generally 
flatter.   The  crossover  point,  where  the  maximum  velocity  at  both  regions 
has  about  the  same  magnitude,  occurs  at  a  flow  rate  of  about  26.5  lbm/sec, 
as  shown  in  Figs.  12,  14,  15,  22,  and  29.   To  this  flow  rate  corresponds 
a  flow  coefficient   Cp  of  about  0.44,  which  is  the  average  of  the  flow 
coefficients  obtained  from  the  calculated  and  measured  axial  velocity 
profiles  presented  in  the  abovementioned  figures.   Except  for  the  regions 
near  the  hub  and  the  tip,  the  axial  velocTty  is  nearly  constant  at  this 
flow  rate.   Further  evidence  that  the  design  criteria  of  a  free-vortex 
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flow  were  most  nearly  met  at  a  (V  of  0.44  is  given  in  Figs.  33  through 
35,  where  it  is  shown  that  total  pressure  is  nearly  constant  except  for 
the  hub  and  tip  regions.  However,  this  uniformity  is  not  so  well  real- 
ized after  the  second  and  third  stage  rotor.  For  this  discussion  it  is 
of  interest  to  note  that  Ref .  4  gives  a  design  value  of  (\)  of  0.45  for 
the  free-vortex  blading  at  rotor  and  stator  stagger  angles  of  43.8  and 
29.8  degrees  respectively. 

Reduced  total  pressures  and  low  axial  velocities  occur  in  the  regions 
near  the  hub  and  near  the  outer  casing  downstream  of  all  rotor  rows. 
They  are  caused  by  end  and  tip  clearance  losses,  by  the  wall  boundary 
layers,  and  by  the  effect  of  the  rotating  hub.   These  conditions  produce 
large  entropy  gradients  at  the  inner  and  outer  radii  of  the  annulus  for 
the  reason  shown  in  the  explanation  of  Eq.  3  in  this  subsection.   Their 
effect  near  the  hub  extends  from  radius  ratios  of  0.6  at  the  hub  to 
about  0.64  for  the  first  rotor,  0.66  for  the  second  rotor,  and  0.69  for 
the  third  rotor.   Near  the  outer  casing,  they  influence  the  velocity 
distributions  at  radius  ratios  of  about  0.91  after  the  first  rotor  and 
about  0.94  after  the  second  and  third  rotor.   The  regions  which  are  af- 
fected by  the  hub  and  tip  losses  become  larger  as  the  fluid  proceeds 
downstream  since  the  losses  accumulate  from  stage  to  stage  and  further 
deteriorate  the  flow  profiles. 

The  velocities  at  the  outer  casing  are  difficult  to  measure.   Since 
the  probes  extended  into  the  flow  annulus  from  the  outer  casing,  radial 
position  changes  of  the  probes  of  the  order  of  0.01  inches  at  the  outer 
radius  caused  changes  in  the  measured  data  of  up  to  40  per  cent.   It  is 
believed  that  these  discrepancies  were  caused  by  small  errors  in  the  probe 
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position  indication.   The  probes  retract  into  the  radial  access  holes, 
and  small  changes  in  radial  position  cause  the  access  holes  to  block 
the  measuring  holes  of  the  probes. 
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SECTION  9 
CONCLUSIONS  AND  RECOMMENDATIONS 

The  accurate  determination  of  the  radial  distribution  of  absolute 
velocity  in  magnitude  and  direction  is  important  for  the  evaluation  of 
the  performance  of  turbomachines.   If  these  distributions  were  known,  the 
velocity  triangles  would  be  completely  determined.   Knowledge  of  the  veloc- 
ity triangles  allows  direct  comparison  between  the  applied  design  criteria 
and  the  actual  results  for  a  turbomachine. 

The  accurate  predictions  of  these  velocity  triangles  is  of  particular 
importance  in  the  design  of  multi-stage  axial-flow  machines.   With  the 
exception  of  the  first  stage,  the  inlet  conditions  of  a  given  stage  are 
the  discharge  conditions  of  the  preceding  stage.   If  a  stage  has  inlet 
conditions  different  from  those  for  which  it  was  designed,  its  perform- 
ance is  likely  to  be  marginal  and  could  produce  discharge  conditions 
that  deviate  even  more  from  the  design  values.   The  result  is  that  the 
machine  will  not  be  able  to  handle  the  desired  enthalpy  change  at  the 
expected  level  of  efficiency.   The  radial  distributions  of  losses  and 
flow  deflections  in  a  row  of  blades  are  a  function  of  inlet  conditions 
as  well  as  the  design  of  the  blades.   A  blade  element  designed  for  a 
certain  range  of  flow  incidence  angles,  but  operating  outside  of  this 
range,  may  either  stall  or  not  develop  the  required  blade  forces,  de- 
pending on  whether  the  incidence  angle  is  too  large  or  too  small.   Optimum 
design  of  multi-stage  axial-flow  machines  therefore  requires  that  the 
discharge  conditions  of  any  stage  match  the  design  inlet  conditions  of 
the  following  stage. 

Since  designs  of  turbomachines  are  based  on  the  results  of  experimental 
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work  and  the  experience  of  the  designer,  checking  the  results  of  new 
designs  is  important  for  two  reasons:   first,  to  evaluate  the  perform- 
ance of  a  new  design,  and  second,  to  improve  the  design  procedures  for 
future  machines. 

The  NISRE  method  represents  a  fairly  accurate  and  simple  method  for 
determining  the  velocity  profile  downstream  of  an  axial-flow  rotor.   The 
data  which  are  required  can  be  obtained  by  small  and  accurate  probes  in 
the  vaneless  space  between  blade  rows  with  small  axial  clearance.   The 
method  as  presented  in  this  study  is  limited  by  the  fact  that  the  ef- 
fect of  meridional  streamline  curvature  is  ignored.   For  designs  where 
the  curvatures  of  the  meridional  streamlines  are  appreciable  ,  the  comput- 
ing process  could  be  improved  in  the  manner  carried  out  by  Harrison  (2). 

Further  experimentation  seems  necessary  to  accurately  predict  the 
gross  weight  flow  blockage  factors,  particularly  with  machines  having 
many  stages.   The  blockage  factors  used  for  the  velocity  profiles  after 
the  first  rotor  were  determined  by  arbitrary  variation  until  accurate 
matching  with  the  measured  profiles  were  obtained.   By  changing  the 
length  of  the  inlet  duct  of  the  O.N.R.  compressor  the  effect  of  boundary 
layer  growth  on  blockage  factors  could  be  investigated. 

The  hot-wire  anemometer  was  not  used  for  an  analysis  of  the  unsteady 
flow  in  the  passages  downstream  of  the  rotors.   The  presented  average 
velocities  were  read  on  an  instrument  which  heavily  damps  the  signal  and 
cannot  follow  the  high  frequency  response  of  the  hot-wire  anemometer 
which  is  up  to  50  kilohertz.   However,  some  oscilloscope  traces  were 
photographed  downstream  of  the  second  rotor.   Sample  pictures  of  these 
oscillograms,  with  the  hot-wire  probe  at  a  radius  of  13  inches,  are  pres- 
ented in  Fig.  78.   The  oscillogram  marked  (a)  was  made  with  the  throttle 
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valve  open,  and  the  oscillogram  marked  (g)  was  made  at  a  throttle  valve 
counter  setting  of  425.   The  measured  flow  rates  are  indicated  in  the 
figures.   The  time  for  successive  blades  to  pass  by  the  hot-wire  is 
1.675  milliseconds,  corresponding  to  a  frequency  of  about  0.6  kilohertz 
which  is  well  within  the  frequency  response  of  the  hot-wire  anemometer. 
Therefore,  the  hot-wire  probe  will  be  able  to  measure  the  flow  condi- 
tions in  the  wakes  between  neighboring  blades. 

The  low  signal  amplitudes  in  Fig.  78,  which  are  marked  by  arrows, 
were  found  to  have  the  time  interval  of  1.675  milliseconds.  Hence,  they 
are  caused  by  blade  wakes.   The  trace  of  the  hot-wire  output  moves  from 
left  to  right  on  the  oscillograms.   Hence,  the  wakes  start  on  the  con- 
cave side  of  the  blade  profiles  and  end  on  the  convex  side  of  the  pro- 
files.  Therefore,  the  flow  profiles  from  left  to  right  between  the 
indicated  wakes  start  from  the  pressure  side  of  one  blade  and  extend  to 
the  suction  side  of  the  neighboring  blade. 

An  interesting  phenomenon  occurs  in  the  oscilloscope  traces  at  a 
flow  rate  of  19.26  lbm/sec  which  is  evident  from  the  oscillogram  marked 
(g) .  An  additional  trace  of  lesser  intensity  but  greater  amplitude  is 
superimposed  on  the  aforementioned  wake  pattern.   This  pattern  is  periodic 
in  nature  and  seems  to  be  an  indication  of  the  propagating  stall  discussed 
by  Iura  and  Rannie  (7). 

If  the  resolution  of  the  oscilloscope  traces  could  be  improved  or  if 
the  hot-wire  signal  were  recorded  on  magnetic  tape  and  transferred  to  a 
strip  chart  recorder,  the  wake  pattern  of  the  flow  profiles  could  be 
analyzed  in  greater  detail  to  detect  incipient  stall  conditions.   Using 
instrumentation  of  this  type,  it  seems  feasible  to  apply  the  momentum 
theorem  to  determine  the  losses  associated  with  either  different  blade 
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shapes  or  various  stagger  angles  of  a  given  set  of  blades.   The  O.N.R. 
compressor  is  suitable  for  such  studies  since  it  is  possible  to  vary 
rotor  stagger  angles  and  to  install  different  types  of  blades.   Such  a 
study  would  require  accurate  calibrations  of  the  hot-wire  probes  to 
enable  measurements  of  velocity  in  absolute  values  rather  than  as  ratios 
with  a  known  velocity.  A  test  rig  is  now  being  designed  at  the  Naval 
Postgraduate  School  for  such  calibrations. 
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FIGURE  3 
COMPRESSOR  INLET  DUCT 


FIGURE  4 
COMPRESSOR    INLET  WITH    [NSTALLED    FLOW    RATE    SURVEY    INSTRUMENTS 


FIGURE    5 
MICROMANOME 
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FIGURE    6 
DETAIL  OF    INLET    HOT-WIRE   ANEMOMETER   PROBE 
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FIGURE  10 
VIEW  OF  COMPRESSOR  SHOWING  FLOW  PROBE  CARRIAGE 
[NSTALLED  AND  ASSOCIATED  MANOMETER  BOARD 
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FIGURE    16 

AXIAL  VELOCITY  RATIO  VS 

RADIUS    RATIO 

FIRST  STAGE   ROTOR 
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FIGURE 
AXIAL  VELOCITY  RATIO 
RADIUS   RATIO 
FIRST   STAGE   ROTOR 
W  =   22.35    lbm/sec 
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FIGURE    18 

AXIAL  VELOCITY  RATIO  VS 

RADIUS   RATIO 

FIRST  STAGE   ROTOR 

W  =   20.75    lbm/sec 
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FIGURE   25 
AXIAL  VELOCITY  RATIO  VS 
RADIUS   RATIO 
SECOND  STAGE   ROTOR 
W  =   20.36    lbm/sec 
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FIGURE   30 

AXIAL  VELOCITY  RATIO  VS. 

RADIUS    RATIO 

THIRD  STAGE   ROTOR 

W  =   24.58    lbm/sec 
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FIGURE   32 

AXIAL  VELOCITY  RATIO   VS. 

RADIUS    RATIO 

THIRD  STAGE  ROTOR 

W  =  20.71  lbm/sec 
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FIGURE   33 
TOTAL  PRESSURE   COEFFICIENT  VS.    RADIUS  RATIO 
FIRST  STAGE   ROTOR 
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FIGURE  3*+ 
TOTAL  PRESSURE  COEFFICIENT  VS.  RADIUS  RATIO 
SECOND  STAGE  ROTOR 
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FIGURE  35 
TOTAL  PRESSURE  COEFFICIENT  VS.  RADIUS  RATIO 
THIRD  STAGE  ROTOR 


76 


R/R0 

FIGURE  36 
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FIGURE   38 

ex'    VS.    RADIUS   RATIO 

FIRST   STAGE  ROTOR 

W  =   26.97    lbm/sec 


FIGURE   37 

cX    VS.    RADIUS   RATIO 

FIRST  STAGE   ROTOR 

W  =  27.90   lbm/sec 
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FIGURE  39 

o<  VS.  RADIUS  RATIO 

FIRST  STAGE  ROTOR 

W  =  26.17  lbm/sec 
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FIGURE  40 
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FIGURE  42 
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cX  VS.  RADIUS  RATIO 

FIRST  STAGE  ROTOR 

W  =  22.35  lbm/sec 
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FIGURE   43 
o<   VS.    RADIUS   RATIO 
SECOND   STAGE   ROTOR 
W  =   28.79    lbm/sec 
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FIGURE  45 
ex1  VS.  RADIUS  RATIO 
SECOND  STAGE  ROTOR 
W  =  2  7.16  lbm/sec 
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FIGURE  44 
oc  VS.    RADIUS   RATIO 
SECOND   STAGE   ROTOR 
W  =   28.09    lbm/sec 
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FIGURE  46 
o<  VS.    RADIUS   RATIO 
SECOND   STAGE  ROTOR 
W  =   26.46    lbm/sec 
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FIGURE  47 
CX  VS.  RADIUS  RATIO 
SECOND  STAGE  ROTOR 
W  =  24.83  lbm/sec 
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FIGURE  49 
oc  VS.    RADIUS   RATIO 
SECOND   STAGE   ROTOR 
W  =   20.36    lbm/sec 
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FIGURE  48 
o<   VS.    RADIUS   RATIO 
SECOND  STAGE   ROTOR 
W  =  22.85    lbm/sec 
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FIGURE    50 

o<    VS.    RADIUS   RATIO 

THIRD   STAGE  ROTOR 

W  =   28.62    lbm/sec 


FIGURE   52 

o<  VS.    RADIUS   RATIO 

THIRD   STAGE   ROTOR 

W  =   27.19    lbm/sec 
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FIGURE    51 

oC   VS.    RADIUS   RATIO 

THIRD   STAGE  ROTOR 

W  =   27.91    lbm/sec 
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FIGURE  53 

oC  VS.  RADIUS  RATIO 

THIRD  STAGE  ROTOR 

W  =  26.08  lbm/sec 


81 


FIGURE  54 

oC  VS.  RADIUS  RATIO 

THIRD  STAGE  ROTOR 

W  =  24.58  lbm/sec 
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FIGURE  56 

o<  VS.  RADIUS  RATIO 

THIRD  STAGE  ROTOR 

W  =  20.71  lbm/sec 
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FIGURE  55 

o<"  VS.  RADIUS  RATIO 

THIRD  STAGE  ROTOR 

W  =  22.05  lbm/sec 
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FIGURE   57 

STATIC   PRESSURE   COEFFICIENT 

VS.   RADIUS   RATIO 

FIRST   STAGE  ROTOR 

W  =  28.49   lbm/sec 


R/R0 


NISRE 

O   FLOW 

Q   HOT-WIRE 


1.0 


0.2 


0.1 


< 

PU 


•0.1 


FIGURE   58 

STATIC   PRESSURE   COEFFICIENT 

VS.    RADIUS   RATIO 

FIRST   STAGE   ROTOR 

W  =   27.90    lbm/sec 
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FIGURE   59 

STATIC   PRESSURE   COEFFICIENT 

VS.    RADIUS  RATIO 

FIRST   STAGE  ROTOR 

W  =   26.97    lbm/sec 
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FIGURE   60 

STATIC   PRESSURE   COEFFICIENT 

VS.    RADIUS  RATIO 

FIRST   STAGE  ROTOR 

W  =  26.17    lbm/sec 
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FIGURE 
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FIGURE   63 

STATIC    PRESSURE   COEFFICIENT 

VS.    RADIUS   RATIO 

FIRST   STAGE   ROTOR 

W   =   20.75    lbm/sec 
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FIGURE   64 
STATIC   PRESSURE   COEFFICIENT 
VS.    RADIUS    RATIO 
SECOND   STAGE   ROTOR 
W  =   28.79   lbm/sec 
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FIGURE   65 
STATIC   PRESSURE   COEFFICIENT 
VS.    RADIUS   RATIO 
SECOND  STAGE   ROTOR 
W  =  28.09   lbm/sec 
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FIGURE  66 
STATIC   PRESSURE   COEFFICIENT 
VS.    RADIUS   RATIO 
SECOND  STAGE   ROTOR 
W  =  2  7.16    lbm/sec 
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FIGURE   67 
STATIC   PRESSURE   COEFFICIENT 
VS.    RADIUS  RATIO 
SECOND   STAGE   ROTOR 
W  =  26.46    lbm/sec 
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FIGURE   68 
STATIC   PRESSURE   COEFFICIENT 
VS.    RADIUS   RATIO 
SECOND   STAGE   ROTOR 
W   =   24.83    lbm/sec 
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FIGURE   69 
STATIC   PRESSURE   COEFFICIENT 
VS.    RADIUS   RATIO 
SECOND   STAGE   ROTOR 
W  =   22.85    lbm/sec 


0.6 


0.7      R/Rq        0.8 


'NISRE     ~" 

O     FLOW 

D    HOT-WIRE 


0.9 


1.0 


89 


0.3 


0.2 


0.1 


PL, 
I 

Pw 


FIGURE    70 
STATIC   PRESSURE   COEFFICIENT 
VS.    RADIUS   RATIO 
SECOND   STAGE   ROTOR 
W  =   20.36    lbm/sec 
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FIGURE    71 

STATIC   PRESSURE   COEFFICIENT 

VS.    RADIUS   RATIO 

THIRD  STAGE   ROTOR 

W   =   28.62    lbm/sec 
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FIGURE    72 

STATIC   PRESSURE   COEFFICIENT 

VS.    RADIUS   RATIO 

THIRD  STAGE   ROTOR 

W  =   27.91    lbm/sec 
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FIGURE   73 

STATIC   PRESSURE   COEFFICIENT 

VS.    RADIUS   RATIO 

THIRD  STAGE  ROTOR 

W  =  27.19   lbm/sec 
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FIGURE    74 

STATIC   PRESSURE   COEFFICIENT 

VS.    RADIUS   RATIO 

THIRD  STAGE  ROTOR 

W  =  26.08   lbm/sec 
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FIGURE   75 
STATIC   PRESSURE   COEFFICIENT 
VS.    RADIUS   RATIO 
THIRD   STAGE   ROTOR  , 
W  =   24.58   lbm/sec 
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FIGURE    76 

STATIC   PRESSURE   COEFFICIENT 

VS.   RADIUS   RATIO 

THIRD   STAGE   ROTOR 

W  =  22.05    lbm/sec 
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(a)    W   =    28.37    Lbm/sec 
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(b)    W    ■     2  7.91    lbm/sec 


PROBE  POSITION 
DOWNSTREAM  OF 
SECOND  STAGE  ROTOR 
RADIUS  -  13  inches 


OSCILLOSCOPE  SETTINGS 
TIME  SCALE    1  millisec/cm 
AMPLITUDE    10  mi  lllvo It/cm 

HOT-WIRE  CURRENT    25  milliamps 


FIGURE  7b 
HOT-WIRE  ANEMOMETER  OSCILLOGRAMS 
SECOND  STAGE  ROTOR 
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(e)  W  -  24.80  lbm/sec 
FIGURE  78  (continued) 
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APPENDIX  A 


DETAILS  OF  FLOW  RATE  CALIBRATION  COMPUTER  PROGRAM 


This  program  uses  the  data  obtained  from  the  hot-wire  anemometer 
surveys  of  the  inlet  to  calculate  volume  flow  rate.   Reference  velocity 
for  the  hot-wire  reference  setting  is  obtained  from  the  inlet  survey 
Pitot-static  probe  mounted  at  the  inlet  centerline.   Inlet  guide  vane 
velocity  is  obtained  from  the  permanently  mounted  Pitot-static  probe 
ahead  of  the  inlet  guide  vanes. 

This  program  assumes  that  the  hot-wire  instrument  is  set  to  90 
per  cent  of  full  scale  for  the  reference  velocity  at  the  inlet  center- 
line  and  that  the  hot-wire  instrument  is  linear. 
Input 
Card  No.       Format        Fortran  Description 

1  12  KKK  Number  of  runs  to  be 

considered. 

2  313,2110  M  Number  of  radial  data 

points  -  symmetric  on 
either  side  of  centerline 
(maximum  99) 
MM  Number  of  readings  of  total 

pressure  and  velocity  head 
taken  at  the  inlet  centerline 
(maximum  99) 
MMM  Number  of  readings  of  total 

pressure  and  velocity  head 
taken  ahead  of  the  inlet  guide 
vanes  (maximum  99) 
99 


IRUN  Run  number 

ISPEED        Rotor  speed  (rpra) 
4F10.4         TI  Temperature  in  the  inlet  (°F) 

TW  Temperature  in  test  cell  (°F) 

BAR  Barometric  pressure  (inches 

Hg.) 
TBAR  Temperature  of  mercury  column 

of  barometer  (°F) 
8F10.4         PT1R  Total  pressure  at  the  inlet 

centerline  (inches  H  0) 
Q1R  Velocity  head  at  the  inlet 

centerline  (inches  HO) 
8F10.4         PT2R  Total  pressure  ahead  of  the 

inlet  guide  vanes  (inches 
H20) 
8F10.4         Q2R  Velocity  head  ahead  of  the 

inlet  guide  vanes  (inches 
H20) 
8F10.4         R  Radii  at  which  data  are  taken 

from  the  centerline  outward  - 
data  points  must  be  symmetric 
on  either  side  of  centerline 
(inches) 
8F10.4         VRL  Values  of  velocity  ratio  of 

hot-wire  anemometer  to  the 
left  side  of  centerline  looking 
downstream  (dimensionless) 


mo 


8F10.4 


VRR 


Values  of  velocity  ratio  of 
hot-wire  anemometer  to  the 
right  side  of  centerline 
looking  downstream 
(dimensionless) 


Output 

The  computer  progran  output  prints  the  input  data  in  the  same  format 
as  they  are  read  in. 

The  calculated  information  is  presented  as  follows: 


Description 

Barometric  pressure 

Inlet  temperature 

Volume  flow  rate 

Mass  flow  rate 

Velocity  at  inlet  centerline 

Velocity  ahead  of  inlet  guide  vanes 

Mass  density  at  inlet  centerline 
..  (volume  flow  rate) 


Units 
(lb/ft2) 

(°R) 

(ft3/sec) 

(lbm/sec) 

(ft/sec) 

(ft/sec) 

(lbm/ft  ) 

(ft3/sec) 

(velocity  ahead  of  inlet  guide  vanes)   (ft/sec) 
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APPENDIX  B 


DETAILS  OF  FLOW  PROPERTIES  COMPUTER  PROGRAM 


Input 

Card  No. 

Format 

Fortran 

1 

12 

KKK 

This  program  uses  the  total  pressure,  total  temperature,  velocity 
head,  velocity  ratio  from  the  hot-wire  anemometer,  and  absolute  flow 
angle  data  obtained  from  the  radial  surveys  downstream  of  rotor  blade 
rows  and  calculates  flow  properties  based  on  the  NISRE  solution,  flow 
probe  velocity  and  hot-wire  anemometer  velocity.  The  NISRE  solution 
also  requires  total  pressure  and  velocity  head  data  from  the  permanent 
Pitot-static  probe  ahead  of  the  inlet  guide  vanes  in. order  to  calculate 
mass  flow  rate.  • 


Description 
Number  of  sets  of  data  to  be 
read. 

Run  number 

Number  of  radial  data  points 
(maximum  50) 

Number  of  radial  points  re- 
quired from  NISRE  solution 
(maximum  100) 

Barometric  pressure  (inches 
Hg.) 

Temperature  of  mercury  column 
of  barometer  (°F) 
Rotor  speed  (rpm) 
Assumed  gross  weight  flow 
blockage  factor  (BKG  ^  1) 


3I3,7F10.3 


IRUN 
MD 


BAR 

TBAR 

RPM 
BKG 
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6F10. 3,313 


GAM 

UGAS 
FARB 


VA1 


ACC 


TW 

TI 
RO 

RI 


Specific  heat  ratio 
(dimensionless) 
Gas  constant  (ft-lb)/lbm-°R) 
Arbitrary  multiplication 
factor  for  reference  axial 
velocity  in  NISRE  solution, 
renuired  to  start  Newton 
Kapha on  iteration 
(1<  FARB<  1.2) 
Reference  axial  velocity  at 
root  mean  square  radius  for 
NISRE  solution,  which  must 
be  estimated  within  50  per 
cent  of  true  value  (ft/sec) 
Accuracy  to  which  the  NISRE 
solution  is  to  match  the 

actual  flow  rate. 

/Kryn          1  Flow  rate  calculated 
(ACC     1-"A— L T~rT T 

Actual  flow  rate 
0<  ACC  £.  .00001) 
Temperature  in  the  test  cell 
<°F) 

Inlet  air  temperature  (°F) 
Outer  radius  of  the  annulus 
(inches) 

Inner  radius  of  the  annulus 
(inches) 
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MMM 


NREF 


8F10.3 


NS 
PTIR 


QIR 


Number  of  readings  of 
total  pressure  and  velocity 
head  taken  ahead  of  the  in- 
let guide  vanes  (maximum  50) 
The  number  of  the  data  point 
at  which  the  hot  wire  anemome- 
ter was  set  to  90  per  cent, 
counting  from  the  hub  outward 
Stage  number 

Total  pressure  ahead  of  the 
inlet  guide  vanes  (inches 
(H20) 

Velocity  head  ahead  of  the 
inlet  guide  vanes  (inches 


H20) 

The  remaining  input  data  are  from  the  radial  surveys  downstream  of  a 
particular  rotor.   Data  must  be  entered  into  the  program  from  the  hub 
outward,  one  card  per  radial  station. 
5  6F10.4  PTR 


QR 

TTR 


Total  pressure  (inches  H  0) 
Velocity  head  (inches  HO) 
Total  temperature  (milli- 
volts determined  from  an 
iron-constantan  thermocouple 
with  a  32°F  reference  junc- 
tion) 


108 


ALPHA         Absolute  flow  angle,  taken 
positive  from  the  axial 
direction  to  the  direction 
of  rotation  (degrees) 

VR  Velocity  ratio  from  hot-wire 

anemometer  (dimensionless) 

R  Radius  (inches) 

Output 

v 
The  computer  program  output  gives  the  input  data  in  the  same  format 

as  read  in  for  the  purpose  of  checking  input  data  for  errors. 

The  major  part  of  the  output  is  divided  into  three  sections,  the 
NISRE,  flow  probe,  and  hot-wire  calculations  respectively.   The  calcu- 
lated data  for  each  of  the  three  methods  are  presented  in  dimensional  and 
in  dimensionless  form.   The  method  used  to  obtain  the  dimensionless 
quantities  is  described  in  Section  8. 
Dimensional  Form 


Name 

RADIUS 

PT 

P 

TT 

T 

V 

VA 

vu 

VMACH 


Description 
(inches) 

Total  pressure  (psfa) 
Static  pressure  (psfa) 
Total  temperature  (°R) 
Static  temperature  (°R) 
Absolute  velocity  (ft/sec) 

Axial  component  of  absolute  velocity  (ft/sec) 
Peripheral  component  of  absolute  velocity  (ft/sec) 
Mach  number  of  absolute  velocity  based  on  calculated 
static  temperature  (dimensionless) 
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W  Relative  velocity  (ft/sec) 

WA  Axial  component  of  relative  velocity  (ft/sec) 

WU  Peripheral  component  of  relative  velocity  (ft/sec) 

WMACH  Mach  number  of  relative  velocity  based  on  calculated 

static  temperature  (dimensionless) 
U  Peripheral  velocity  (ft/sec) 

ALPHA  Angle  of  absolute  velocity  -  positive  from  the 

axial  direction  to  the  direction  of  rotation  (degrees) 
BETA  Angle  of  relative  velocity  -  positive  from  the 

axial  direction  to  the  direction  of  rotation  (degrees) 
Dimensionless  Form 

All  velocities  are  made  dimensionless  by  dividing  by  the  peripheral  tip 
velocity  (U  ) .  All  temperatures  are  made  dimensionless  by  dividing  by 
the  compressor  inlet  temperature. 


R/RO 


DP/RHO*UT  **2 
PHIAV 


Radius  ratio  -  Radius  divided  by  tip  radius 


DPT/RHO*UT  **2      (P  -P  )/(%PU_  ) 

t  a    v   1 

2 


(p-pfl)/(%  pu,p 

Numerical  average  of  the  dimensionless  axial 
velocity  components. 
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